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Abstract 

The phenomenon of the finite-temperature induced quantum numbers in fermionic 
systems with topological defects is analyzed. We consider an ideal gas of twodimen- 
sional relativistic massive electrons in the background of a defect in the form of a 
pointlike magnetic vortex with arbitrary flux. This system is found to acquire, in addi- 
tion to fermion number, also orbital angular momentum, spin, and induced magnetic 
flux, and we determine the functional dependence of the appropriate thermal averages 
and correlations on the temperature, the vortex flux, and the continuous parameter of 
the boundary condition at the location of the defect. We find that nonnegativeness of 
thermal quadratic fluctuations imposes a restriction on the admissible range of values 
of the boundary parameter. The long-standing problem of the adequate definition of 
total angular momentum for the system considered is resolved. 



1 Introduction 

Quantum fermionic systems in different nontrivial topological backgrounds (kinks, vortices, 
monopoles, skyrmeons etc.) can possess rather unusual properties (e. g. fractionization of 
quantum numbers) [HE], see reviews in Refs.jSlllj, and an interest to finite-temperature 
effects in such systems jS] has been recently revived [HI El E] • in particular, planar sys- 
tems with a topological defect in the form of a pointlike magnetic vortex deserve a thorough 
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examination, since they may be relevant for the description of some condensed matter phe- 
nomena, including superfluidity and superconductivity (TU1 HH C2] , as well as have various 
applications in particle physics, cosmology, and astrophysics |131 IT4] . On the other hand, 
these systems can be of a certain conceptual importance, providing a field-theoretical mani- 
festation of the famous Bohm-Aharonov effect ^3]: they involve second-quantized fermions 
interacting with a vector potential which is caused by a magnetic flux from the inaccessible 
for the fermions region. 

A study of quantum numbers which are induced in the Bohm-Aharonov manner (i. e. by 
a vector potential of a magnetic vortex) started in Refs. (THllIZj- It was shown for a particular 
choice of the boundary condition at the location of the defect that electric charge JH] , mag- 
netic flux and angular momentum [20] are induced in the vacuum of quantized massive 
fermions. The induced vacuum quantum numbers under the most general set of boundary 
conditions which are compatible with self-adjointness of the pertinent Dirac Hamiltonian 
were obtained in Refs. |211 1221 I23| . The finite-temperature induced charge was examined in 
Ref . |24j . Following this line, we consider other finite-temperature induced quantum numbers 
in the present paper. 

We start with the operator of the second-quantized fermion field in a static background, 
*(x,f)=^ e-^(x|A)a A + ^ e -"^<x|A>&+, (1.1) 

(E x >0) (E x <0) 

where a\ and a\ (b\ and b\) are the fermion (antifermion) creation and destruction operators 
satisfying anticommutation relations, 

[a A ,a+] + = [6 Al 6+] + =(A|A , >, (1.2) 
and (x|A) is the solution to the stationary Dirac equation, 

F(x|A) = £ A (x|A) , (1.3) 
H is the Dirac Hamiltonian, A is the set of parameters (quantum numbers) specifying a 

.. — — , t — — 

discrete and the integration (with a certain measure) over continuous values of A. Ground 
state |vac> of the second-quantized theory is defined as 

a A |vac>= &A|vac>= . (1.4) 

Let J be an operator commuting with the Hamiltonian in the first-quantized theory, 

[J,#]_ = 0. (1-5) 

In the case of unbounded operators, commutation of their resolvents is implied, or, to be 
more specific, it is sufficient to require that operators H and J have a common set of 
eigenfunctions, i.e. relation 

J < x|A >= j\ < x|A > (1.6) 
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holds as well as Eq. (J1.3|) . Eigenfunctions < x|A > satisfy the conditions of completeness 
and orthonormality; in general, normalization to a delta function is implied. Thus, in the 
second-quantized theory, the operators of the dynamical variables (physical observables) 
corresponding to H and J can be diagonalized: 



P° = I / d d x\[S>+(x,t), fifc*(x,t)]_- [^ + (x,t), *(x,t)]_ 



;a A - b+b x - -sgn(£' > 



and 



M= l -jd d x [¥+(x,t), J*(x,t)]_ = ^ 



J A 



a A fl A - &a & A - -Sgn(E A ) 



;i-7) 



(1.8) 



d is the space dimension. 

Let us define partition function 



Z(j3,iu) = Spexp -(3(P -fijM) , /3 = {k B T)~ l 



;i.9) 



where T is the equilibrium temperature, ks is the Boltzmann constant, fij is the generalized 
chemical potential, and Sp is the trace or the sum over the expectation values in the Fock 
state basis created by operators in Eq. (jl.2|) . Although this sum becomes divergent in the 
limit of infinite space volume, this will not bother us, since the partition function plays 
a merely supplementary role. The quantities of physical interest are obtained by taking 
derivatives of In Z(/3, /ij), and these latter may appear to be finite in the infinite volume 
limit. 

In particular, one can define an average of operator M over the grand canonical ensemble 



m) 



1 d 



(3dfMj 
Computing averages 



\nZ(P,fij) = Z- l {P,f,j)SpMexp -(3 (P° - fxjM 



{exp^-zij^ + l}- 1 , E x >0 
{exp[/3(-£ A + fi jJx )} + l}- 1 , E\ < 



;i.io) 



;i.n) 



and using the explicit form of P° and M in terms of the creation and destruction operators, 
one gets the spectral integral representation of average (jl.lO|) (see, e.g. Ref.jH]), 



m) 



dErj(E) tanh 



1.12) 



where the appropriate spectral density is 



1 



Tj (E) = -Im Tr J(H — E — iO) 



-i 



;i.i3) 
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Tr is the trace of an integro-differential operator in the functional space: TrU = J d d x tr 
< x|Z7|x >; tr denotes the trace over spinor indices only; note that the functional trace 
should be regularized and renormalized by subtraction, if necessary. 

Taking J — I, where I is the unit matrix in the space of Dirac matrices, one gets 
M = N, where N is the fermion number operator in the second-quantized theory, then /i/ is 
the usual chemical potential. In the d — 1 case fermion number is the only observable which 
is conserved in addition to energy. In more than one dimensions there are more conserved 
observables. In particular, in the d = 2 case, in addition to energy and fermion number, also 
total angular momentum is conserved when the system is rotationally invariant. 

Now let us consider an observable which is not conserved and denote the appropriate 
operator in the first-quantized theory by Q. Then the corresponding operator in the second- 
quantized theory, 



O 



;i.!4) 



is not diagonalizable. Nevertheless, its average over the grand canonical ensemble can be 
defined in a manner similar to Eq. ()1.10|) . 



6) =Z^((3,fij)Spdexp -f3(P -^M 
One can get an appropriate spectral integral representation, 



;i.i5) 



o 



dEr^i(E) tanh 



;i.i6) 



where 



tq(E) = -ImTrVtiH — E — iO) 

7T 



-1 



;i.i7) 



In the present paper we shall be dealing with the averages over the canonical ensemble: 



M(T) 



M 



0{T) = (6 



'1.181 



In addition to them we shall be considering also such quantities as the correlation of the 
conserved and nonconserved observables 



A(T;0,M) = (6m) -(d) (m) 

and the quadratic fluctuation of the conserved observable 



1.19) 



A(T; M, M) = (M 



m) 



;i.2o) 



Using Eqs. p,12|l and (jl.lfij) . one can get the spectral integral representation for Eqs. p,19j 
and (IT^ : 



A(T;6,M) 



d_ 
(3 \d/jj 



O 



l - I dE tqj(E) sech 2 (^E 



1.211 
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and 

oo 
— oo 

where the appropriate spectral densities are obtained from Eqs. f)1.17j) and (jl.l3|) by inserting 
an additional power of operator J into the trace. 

2 Observables of the planar fermionic system 
in the background of a magnetic vortex defect 

We consider a spinor field which is quantized in the background of a static magnetic field in 
2 + 1-dimensional space-time. The Dirac Hamiltonian takes form 

H = -ia[d-ieV(x)]+/3m, (2.1) 

where V(x) is the vector potential of the field strength -B(x) = d x V(x). The Clifford 
algebra in this case has two inequivalent irreducible representations which can be differed in 
the following way: 

a l a 2 (3 = is, s = ±1 . (2.2) 
Choosing the j3 matrix in the diagonal form, 

P = 0-3 , (2.3) 

one gets 

a 1 = -e^ 3Xs a 2 e~^ 3Xa , a 2 = se^ 3Xs a x e~i a3Xs , (2.4) 

where a x , a 2 , and u 3 are the Pauli matrices, and Xi an d X-i are the parameters varying in 
interval < Xs < 27r to go over to the equivalent representation. Note also that in odd- 
dimensional space-time the m parameter in Eq. (|2.1|) can take both positive and negative 
values; a change of sign of m corresponds to going over to the inequivalent representation. 

If a magnetic field is invariant under rotations of the twodimensional space around its 
origin, then one has 

(x x 3) [d x V(x)] = , (2.5) 
and a generator of rotations takes form 



r 

J = -jxx [d - z'eV(x)] + ^sp + e J drr[d x V(x)] 



(2.6) 



r = (x 1 ) 2 + (x 2 ) 2 and tp = arctan(x 2 /x 1 ) are the polar coordinates. One can easily verify 
that operator J (j2.6|) commutes with operator i/ (j2.1J) . 

In the right hand side of Eq. (J2.6j) . the first two terms represent the orbital and the spin 
parts of the angular momentum of the charged matter field, whereas the last term represents 
the angular momentum of the background field. In the nonsingular long-range gauge 

x-V(x)=0, (2.7) 
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one gets 

x x V(x) = I drr[dx V(x)] , (2. 



and Eq. (J2.6j) takes form 

J = -ixxd+^sp. (2.9) 

The above is relevant for the case of an extensive configuration of the background mag- 
netic field (see, e.g., Ref.j2SJ). Turning now to the case of the background in the form of a 
magnetic vortex defect, let the central region (e.g. a disc of radius 5) be impenetrable for 
the charged matter and the background field strength be nonvanishing only in this region 
(i.e. the region of the defect). Then the angular momentum operator outside the central 
region consists of two parts, orbital and spin, 

J = -ixx [0 - ieV(x)] + . (2.10) 

As is well known (see, e.g., Ref.[15j), due to nonvanishing flux of the background field in the 
inner region, 



$ = J drr[dx V(x)] , (2.11) 



the vector potential cannot be made vanishing everywhere in the outer region. In particular, 
in the gauge (|2.7|) one gets at r > 5: 

V 1 ^) = -Qr^smp, V 2 (x) = ^r' 1 cosp , (2.12) 

and Eq. (j2.1(J|) takes form 

J=-ixx9-e$+^. (2.13) 

Thus, contrary to the case of the extensive background field configuration when the 
angular momentum is quantized in half-integer values, 

3= n +\> ( 2 - 14 ) 

(this is evident in the gauge (|2.7)L see Eq. (|2.9j) ). in the case of the vortex defect, the angular 
momentum is quantized in units 

j = n+~-e$, (2.15) 

and this results in such fascinating quantum-mechanical concepts as anyons and fractional 
statistics [2B]. However, in the latter case one can take operator 

J' = _ix x [d-zeV(x)] + ^sp + E, (2.16) 



6 



as well as an operator of conserved quantity; here H is an arbitrary constant. In partic- 
ular, choosing S = e<3>, one gets in the gauge (|2.7|) the same expression as Eq. (j2.9|) and, 
consequently, half-integer eigenvalues. The arguments in favour of such a definition of the 
angular momentum operator are given in Ref.j2ZJ. Not going into details of the discussion 
at the quantum-mechanical level, we would like to emphasize that the problem of the proper 
definition of the angular momentum operator might be resolved in the framework of the 
second-quantized theory at nonzero temperature. Indeed, the quadratic fluctuation (jl.22|) 
of the physically meaningful observable has to be nonnegative, and this, as we shall see in 
Section 6, imposes a definite restriction on the choice of the appropriate operator in the 
first-quantized theory. 

Turning to the nonconserved observables, it is natural to consider, in the capacity of Q, 
the orbital angular momentum operator, 

A = -ixx [d - zeV(x)] , (2.17) 

and the spin operator, 

£ = ^/3. (2.18) 
In addition to these, we consider also operator 1 

e 2 

VL = — xxck, (2.19) 

4:71 

which corresponds to the observable with the physical meaning of the induced magnetic flux 
multiplied by e. Really, the latter quantity at finite temperature is 

oo 

0(T) = e J drrB {I \r) , (2.20) 
'<) 

where the induced magnetic field strength B^ is rotationally invariant and satisfies Maxwell 
equation 

d r [rB {I) (r)] = xxj(x) , (2.21) 

with induced current 



i(x) = — tr / x - , 

JV ' 2 \ V 2 

Solving Eq. (|2.21|) and substituting the solution into Eq. (j2.20|) . one gets 



a tanh I —f3H 



(2.22) 



0(T) = £ J rf 2 x[xxj(x)], (2.23) 

r>8 

which is the thermal average of operator Q (J1.14|) constructed from ^ ()2.19|) . 

In the present paper we shall compute thermal averages of the above observables, as 
well as correlations of conserved and nonconserved observables and quadratic fluctuations of 

1 Note that e 2 has dimension of mass in the case of 2 + 1-dimensional space-time. Hence O (|2 . 1 Q|) is 
dimensionless, as well as A (|2.17|) and S (I2.18|) . 
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conserved observables; note that the average and fluctuation of charge (i. e. fermion number 
times e) have been computed earlier [21] . Thermal characteristics of the quantized fermionic 
matter in the background of a magnetic vortex defect depend both on vortex flux (j2.11|) and 
the parameter of the boundary condition for the matter field at the edge of the defect. 



3 Spectral densities and traces of resolvents 

In order to compute thermal characteristics one has to determine spectral densities t^{E), 
tz(E), Ta(E), t-zj(E), Tj2(E), tqj(E), which are imaginary parts of the appropriate func- 
tional traces, see, e.g., Eqs. (jl.l3J) and (jl.l7|) . Actually, integrals over the real energy spec- 
trum can be transformed into integrals over a contour on the complex energy plane, thus 
yielding a representation of thermal characteristics through the traces directly. In particular, 
one gets 

M(T) = -~ J tanh (~/3uA Tr J{H - u)' 1 , (3.1) 
c 

and 

A(T; M, M) = 1 J S ech 2 Q/fc,) Tr J 2 (H - , (3.2) 

c 

and similarly for other characteristics; here C is the contour (— oo -HO, +oo -HO) and (+oo — 
iO, — oo — iO) in the complex cu-plane. 

The kernel of the resolvent (the Green's function) of the Dirac Hamiltonian in the coor- 
dinate representation is defined as 

G>, <p- r', <p') = (r, <p\(H - u)-y, <p') . (3.3) 

Using Eqs. ([2.3|) and (|2.4j) . one can expand Eq. (j3.3|) in modes in the following form 

G»(r in- r' to') - — V e m ^-^ ( ^ ^ ^ r')^'"^ \ ( o 4) 

G [r, <p, r , if ) - ^ ^ e ^ ^ r ,y [s ^ Xs) ^ ^fc-V) ) • ^.4j 

In the background of magnetic vortex defect (|2.12j) . Hamiltonian (|2.1|) takes form 

H = -ia r d r - ir~ 1 a' p (d v - ze$) + (3m , (3.5) 

where 

a r = a 1 cos ip + a 2 sin ip , d 9 = — a 1 simp + a 2 cos<p . (3.6) 

If a size of the defect is neglected (8 — > 0), then a parameter of the boundary condition at 
the location of the defect (at r = 0) exhibits itself as a parameter of a self-adjoint extension 
of the Hamiltonian operator. Partial Hamiltonians are essentially self-adjoint for all n, with 
the exception of n = no, where 

no = [e$] + (3-7) 
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{uj is the integer part of quantity u (i.e., the largest integer which is less than or equal 
to u). The partial Hamiltonian for n = n requires a self-adjoint extension according to 
the Weyl-von Neumann theory of self-adjoint operators (see, e.g., Ref.j2H])- Appropriately, 
radial components a n , b n , c n , and d n in Eq. (j3.4j) with n ^ n are regular at r — > and 
r' — > 0, whereas those with n = n Q satisfy conditions (for details see Ref. [21]): 

cos + h m (|m|r) F a no (r; r — — sgn(m) sin (^ s ~^ + lhn(|m|r) 1 ~' F 5 no (r; r ') 
cos + hm(|m|r) F (i no (r; r') = — sgn(m) sin (^ s ~^ + ^ lim(|m|r) 



and 

cos ( s— H — J lim(|m|r') F a no (r; r') = — sgn(m) sin ( s— H — J lim (Imlr') 1- ^^^; r') 
\ 2 4/ r'— >o " \ 2 A J r'— »0 

cos + ^ lim(|m|r') F 6 no (r; r') = — sgn(m) sin (^ s ~^ + ^ lim(|m|r') 1_F c no (r; r') 

(3.9) 

where 

'1, u > 
-1 , it < 

is the self-adjoint extension parameter, and 



sgn(w) 



F = s{m + I _ I s > (3.10) 

{[m]} = u — [uj is the fractional part of quantity u, < {[u]} < 1; note here that Eqs. (j3.8|) and 

(13. 9j) imply that < F < 1, since in the case of F = s all radial components obey the 

condition of regularity at r — > and r' — > 0. Note also that Eqs. ()3.8|l and ()3.9|) are periodic 
in with period 27r. 

The radial components of the resolvent kernel have been determined in Ref . |24| , and we 
list them in Appendix A. 
Let us consider quantities 

oo 

/oo 
diptr [AG u (r^;r'^)] = [(n - e$)a n (r; r') + (n + s - e$)c n (r; r')] , (3.11) 

M n=— oo 

„ oo 

y ^tr[£G>,^r>)] = - S ^ [a n (r; r0 - ^(r; r')] , (3.12) 



o 



^tr[fiG w (r,^r',^)] = — sr ^ [6 n (r; r') + d n (r; r')] , (3.13) 

n n=— oo 

oo 

J dvtr[Y,JG"{r ) ^r\ V )] = l -s ^ fn - e$ + [a n (r; r') - c n (r; r')] , (3.14) 

n — — fYi > ' 



7 00 / i \ 2 

/ cbptr[j 2 G u (r,<py,<p)]= U-e$ + -aJ [a n (r; + ^(r; 0] , (3-15) 

q n=— oo ^ ' 

G^tr[a/<^(r,^;r>)] = |-sr ^ f n - e$ + -sj [6 n (r; r') + d n (r; r')] , (3.16) 

where operators J, A, E, are given by Eqs. (j2.13j) . (|2.17|l - (|2.19jl . correspondingly. Using 
the explicit form of a n , b n , c n , d n given in Appendix A, we perform summation over n and 
get in the case of r' > r and Im k > \Re k\ (see Appendix B): 

2lT oo 

/ r . ™ ,/ / m ssin(F7r) f, ( n 2 rr' r 2 + r' , s T ^ „ S1 

d^tr [AG u (^;r»] = ^^J d V exp (~^ 2^ V ) i K ^y)- K ^Mv)}+ 

o o 
2 s sinf Ftt) 

+— i >\-F(uj + m) t&n^K F Ur)K F Ur') + (l~F)(u;-m)e lFn K 1 ^ F (Kr)K 1 _ F (Kr' 

^(tan^ + e li,7T ) L 



(3.17) 



2tt °° 2 2 2 

dptr [S G"^(r, y>; r', <p)] = sm K (n\r-r'\)- m J d V ex P — V 



oo 

3 y 



x ; £_ / rf Me - w [(l - F)K F (u) + F^_ F (n)] - (2F - l)[^(y) - V + 



u 



+ ssmCFvr) , + ^ ^t^^t^ _ t u _ m yF, Kl _ F ( Kr \ Kl _ F < Kr >\\ 
7r(tani/ w + e^ 71 ") L J 

(3.18) 



2tt 

/r^^v, ,/ / \ i e 2 ssm(FTt) f, ( K 2 rr' r 2 + r' \ rr ^ . . Tr , 
d^r[OG"(r,^r>)] = ^T^j^ eXp (~^ 2^7 [ 

o o 

+ o 27r^T^ Kr [taa^^i-FW^K) - e^H^i-H^O] , (3.19) 
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2-k 



jd^tr [EJG^(r,^r\cp)]^K (K\r-r'\)- S -^^ J dy exp 



K 2 rr' r 2 + r' 2 



2rr 



x 



4F(1 -F)y % 
y 

(F - |) sin(FTr) 
7r(tanz/ w + e ii?7r ) 



e~ u [{l-F)K F {u) +FK 1 _ F (u)} - 



(F - I) . 



+ m 



2y 



— y x 



2F(1 - F) 

[(u+m) tanv u K F (Kr)K F (Kr') - (uo - m)e lF7T K 1 _ F (Kr)Ki_ F (Kr')] , 



(3.20) 



2n 



J d^tr [j 2 G>,^;r»] =| 



Krr 



Kq(k\t — r'\) + 4, Ki(/t|r — r'l) 

r — r 



sin(F 



7T 



7T 



o?t/ exp 



k tt r + r 



2y 2rr 



— y x 



X 



2F(1 



4 +y y F(l-F) 



+ m 



- ^^(y)] - - F)Kp(y) + FX 1 _ F (y)] ^ + 



2(F - -)sinfF7r) 

+ 2; . [(w + m) tanz/ w ir F (/tr)ir F (/tr / ) + (w - m)e^Ki_ F (/tr)fs: 1 _ F (/tr / )] , 

(3.21) 



7r(tan + e iFw ) 



2tt 



/e — | — 7^ ' 
diftr [f2JG w (r, ip; r', </?)] = — nr- KAn\r — r'l) — 

47T |r — r'\ 

o 

e 2 sin(F7r) 



K 2 rr' r 2 — r' 2 
dy l + ^^4 



K 2 rr' r 2 + r' 2 



— y x 



8F(1 - F)tt 2 J V V 4rr ' 7 V 2 ?/ 2rT 
o 

x \ e yj due- u [{l-F)K F {u) + FK l . F {u)}-{2F-l)y[K F {y)-K l _ F {y)} U 



+ 



e 2 sin(F7r) 
4^2 



rfy exp 



K 2 rr' r 2 + r /2 



2y 2rr 



^y) [(l-Fjif^t/j + F^)]- 



e 2 (F - \) sin(FTr) 
27r 2 (tanz/ a; + e iFv ) 



nr [tanz/ CJJ f^ 1 _ F (/^r)i^ F (/^r , ) - e^ir F (/tr)K 1 _ F (/tr')] , (3.22) 
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where k = —ik, K p (u) is the Mackdonald function of order p, and tanz/ w is given by 
Eq.flXTTHJ). The first terms in the right hand sides of Eas. (|3~T5j) . (pT2njl - (pT22|l diverge in 
the limit r' —> r. These terms coincide with expressions corresponding to the case of absence 
of the vortex defect (see also Appendix B): 



2tt 



diptr [S G w (r, ip; r', cp)) | e $ =0 = smK (K\r — r'|) 



(3.23) 



2tt 



J dtp tr [E J G u (r, <p; r', if)} | e4=0 = |# («k -r'\) 



(3.24) 



2tt 



/ diptr [j 2 G u {r,ip;r',(p)\ | e$=0 = ^ 



K («k - r 'l) + 4 1 jr K x {k\t - r'\) 



(3.25) 



2?r 



/c r I r 1 
diptr [QJ G w (r,^;r',y?)] U=o = ^ Kr i 77^i( K l r - r '\) 
An r — r 



(3.26) 



and, otherwise, 



2tt 



c^tr[AG>,y?;r>)]| e$= o = / <fy>*r [fi G>, y>; r', p)] | e$=0 =0. (3.27) 



Thus, quantities (j3.17|) - ()3.22|) are made finite in the limit r' — > r by subtracting expressions 
corresponding to the case of absence of the vortex defect. Integrating over radial variables, 
we get the renormalized traces: 



oo 2tt 



TrK(H-uo)- l = J drr J d V tr [AG u (r^;r^)] 



o o 



lu 2 — m 2 



F 2 (u + m) tanz/^ — (1 — F) 2 {uj — m)e 
tan u u + e iFir 



iFn 



§(F-l|F(i-f> 



(3.28) 



OO 2-7T 



TrT,{H-uj)- l = J drr J dip {tr [S G>, <p; r, ip)\ - tr [E G"(r, <p; r, <p)\ | e0=o } 



o o 



1 s 

2 tu 2 — m 2 



F(c<j + m) tan z/ w — (1 — F)(o; — m)e 
tan z/. + e ii?7r 



F(l - F)m 



(3.29) 
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OO 2-7T 



TrVt{H -u)- 1 = J drr J dcptr [Q G w (r, <p; r, <p)] 



o o 



e 2 sF(l-F) (1 + F) tan^-(2-F)e^ 



67r uj 2 — m 2 



tan i/^ + e 



TrEJ(H-u 



2 w 2 - m 2 



OO 27T 

drr J dip {tr [ZJ G w (r, ip;r, ip)} - tr [£ JG w (r, ^; r, y>)] U =0 } = 
o o 



1 F-i F(a; + m)tanz/ w - (1 - F)(uj - m)e lF7T 1 F(l — F) 



tan i/ w + e iF?r 



+ - 



A uj 2 — m 2 



(3.31) 



Tr J 2 (H -uj)- 1 



2tt 



drr J dp {tr [J 2 G? (r, (p',r, <p)] - tr [J 2 G"(r, p; r, (p)] | e * =0 } 



(F - |) 2 F(u; + to) tan v u + (1 - F) (uj - m)e iFn 



to 2 — m 2 



tan v w + e iFw 



+ 



1F(1-F) 
2 u 2 - m 2 



1 



-F(l-F) 



+ 



w _l(i?-±)ro^ (3.32) 



oo 2tt 



TrttJ(H -co)- 1 = J drr J dip {tr [QJ G w (r, ip;r, <p)] - tr [VtJG u {r, p; r, p)] | e$=0 } 



o o 



e 2 F(l - F) F(l + F) tanz/ w + (1 - F)(2 - F)e 



?7r cj 2 — m 2 



tan z/. + e iF?r 



(3.33) 



where the integration is performed at Ren > \Iuik\ and, then, is continued analytically to 
half-plane Ren > (Irak > 0) which corresponds to the whole plane of complex uj. For 
completeness, we present here the result of Ref. [24 : 

OO 2-7T 

Tr(H-uj)- l = j * drr J dip [tr G w (r, <p\ r, (p) - tr G^(r, p\ r, (p) | e $ =0 ] = 
o o 



uj 2 — m 2 



F{uj + m) tanz/^ + (1 — F)(u — m)e 
tan 24, + e iF7r 



iFn 



-F(l-F)u 



. (3.34) 



There are remarkable relations among different traces. In particular, summing Eqs. (J3.28|) 
and (EH we get 



Tr J(H - uj)- l — —s [ F — - \ Tr (H — uj)' 1 + 



, sF(l-F) 



2 J uj 2 — m 2 



1 1\ 1 

- IF \ oj H — m 

3 V 2/ 2 



(3.35) 
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Trace (|3.3U|) is expressed through traces (|3.29|) and (|3.34|) : 



^ Tr n(H - uj)- 1 = - Tr (H - uj)- 1 - ( F - - ) Tr E{H - u)- x + 



4 



We list also some other relations: 

2 



+ 



1 sF(l-F) 
4 uj 2 — m 2 



(3.36) 



TtJ\H-u)- 1 =[F~ 1 -\tt{H-uj)- 1 + 1 - F{ } F } 

2 uj z — m z 



Tr \J(H - uj)' 1 — —s [ F ) Tr A(H - oj)" 1 - 



F(l - F)u> - ~ (f - \ ) m 



[1 — F(l — F)]u; -f [F--)m \ , (3.38) 



(3.37) 



TrZJ(H-u)- 1 = -s(f- TrECH-u;)-^ 



i, 1^(1-^) 



^2 _ m 2 



, (3.39) 



Trnj(H-uj)- 1 = -s[F--) Tr Sl(H - a;)" 1 + ^— F<K } F } 

1 2 v 7 Y2n uj 2 -m 2 



1 + -F(1-F) 



. (3.40) 



4 Averages 

Similar to Eq. (j3.1j) . we get the thermal averages of orbital angular momentum 

L(T) = ~\j ^~ tanh Tr A(H - uj)-' , (4.1) 

c 

and spin 

(4.2) 



c 



where C is the contour (— oo + ?0, +oo + i0) and (+oo — zO, — oo — iO) in the complex tu-plane. 
Using Eqs. (j3.28|) and (J3.29|) and deforming the contour around the cuts and poles on the 
real axis, we obtain the following expressions for the averages as real integrals: 



/m . sin(F7r) f du , f 1 „ , > 

L(T) = s — / — tanh -dmVu + 1 x 

o 



F 2 u F A + (1 - F) 2 u l - F A- 1 + u { [\ - F(l - F)] (u F A + u^A' 1 ) - (2F - 1) cos(Ftt)} 

[u F A - u l - F A- 1 + 2 cos(Ftt)] 2 + 4(u + 1) sin 2 (F7r) + 

s r , ^ n fl - 2F(1 - F)\E BS + (2F - l)m , /l „ \ . . 
+ i [l-s gn( ^)]l l (2f _^£ + l m tanh(-/3^j , (4.3) 
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and 

oo 

„ frrf . sin(i<V) f du , (\ „ , \ 

S(T) = -s — i '- / — tanh -fimV^Tl x 

27T J uvu+1 \2 J 

o 

Fu F A + (1 - F)u 1 - F A- 1 + u [\u F A + \u 1 ~ F A' 1 - (2F - 1) cos(Ftt)] 
X \u F A - u l - F A- 1 + 2 cos(Ftt)] 2 + 4(u + 1) sin 2 (F7r) 

- S 11 - *-<*" pF-l ( )gJ+m G"*») + i^ 1 " F) ta " h (H • < 4 ' 4) 
where 

, 2F r(i-F) / e tt\ 

r(w) is the Euler gamma function, Ebs is the energy of the bound state in the one-particle 
spectrum, which is determined as a real root of algebraic equation (for details see Ref. |22|) 

( l +m -^L ^-i; («) 

note that the bound state exists at cos© < (A < 0), and its energy is zero at A = —1, 
and, otherwise, one has < \Ebs\ < \ m \ an d 

sgn(E BS ) = ~ sgn(m) [ sgn(l + A' 1 ) - sgn(l + A)] . (4.7) 

Summing Eqs. (j4.3|) and (|4.4j) we get the thermal average of total angular momentum, 
which can be written in the form (compare with Eq. (|3.35jl ): 

M(T) = -s(f- ^ J N(T) + |F(1 - F) tanh [^Prnj , (4.8) 

where the thermal average of fermion number (i. e. electric charge divided by e) is given by 
expression (see Ref.|2^j): 



* f—j—, \ sin(F7r) f du 
N(T) = - ' 



tanh ( -3m\Ju + 1 ] 

u + 1 V2 J 



vr J uy/u + 



Fu F A - (1 - F)u 1 - F A- 1 +u[(F- |) (u F A + m 1 "^- 1 ) - cos(Ftt)] 



[u F A - u l ~ F A- 1 + 2 cos(Ftt)] 2 + 4(u + 1) sin 2 (F7r) 



i [l-sgn(A)] tanh Q/3£ flS ) (1.9) 



oo \ 

Note that both L(T) and 5*(T) are infinite ^divergent as integral J — J at half- integer 



values of e$, unless A = or A 1 = 0. However, this divergence cancels in the sum, and 
one gets 

— tanh -/ 
16 V2' 



M(T) | F= i = — tanh ( -(3m ) . (4.10) 
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In the cases of A = and A 1 = expressions (|4.3|) and (J4.4j) simplify 



and 



L(T) = i( F -l±l) tanh f^ru) , 8 = ±s| (mod 2vr) 



5(T) = -| ( F - \ ± ^ J tanh Q/3mj , 6 = (mod27r) 



(4.11) 



(4.12) 



In the limit T — > (/3 — >■ oo) we get the results of Ref . |2~3] : 



L(0) 



s sgn(mj 
2^ 



^ F 2 i;" 1+F A + 1 - 2F(1 — F) + (1 — Ffv- F A~ l 



and 



S(0) 



4 



s sgn(mj 



F(l-F) - 



+ 2 + t; 1 "^- 



dv Fv~ 1+F A + 1 + (1 - F)v~ F A- 1 



(4.13) 



7r y v^-i 
i 



+ 2 + t; 1 -^- 1 

In the high-temperature limit the averages tend to zero 
L(T -f oo) = 



(4.14) 



, , sin(FTr) T(l - F) 1 - 2F(1 - F) 
s sgn(m) ^rr. — : — tan ( s 



s sgn(mj 

and 
S(T -> oo) 



2tt T(l + F) 1-2F 
sin(F7r) T(F) 1 - 2F(1 - F) 
"27T r(2 - F) 2F- 1 










~2 ' 











( \ m \ 


2 H 


■5) 


\k B T 



1-2F 



2F-1 



< F < 



- < F < 1 
2 



(4.15) 



, ,sin(F7r)r(l-F) 1 / tt 

-ssgn(m) ; tan s 1 — 

& K J 4tt T(l + F) 1-2F V 2 4 

sin(FTr) T(F) 1 



, 7T 

-ssgnfm) ) '- — — '—— cot s 1 — 

4tt r(2-F)2F-l V 2 4 



I I \ !- 2F 1 
771 \ 1 

m • 0<F< 2 

771 \ 1 „ 



(4.16) 



In conclusion of this section, let us consider the thermal average of induced flux (times e) 

0(T) = -\ [ ^tanh([3uj)Trn(H - uo)- 1 , (4.17) 
2 J 2ni 



c 

where Q is given by Eq. (j2.19|) . Using trace identity ()3.36|) . we get 

„2 



0(T) 



2ixm 



4 



N(T) 



F - X - J S-(T) + A f f - i J F(l - F) tanh Q/3m 



(4.18) 
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or 



0(T) 



' sF(l - F) J sin(F7F 



4nm 



n J a ■ 

o 

F A _ „,1--F 4-1 



tanh ( -3m\Ju + 1 ) x 

At + 1 V 2 / 



u A - u A - 2mcos(Ftt) 



+ 



+ - [1 - sgn(A)] 



- m 1 -^^ 1 + 2 cos(Fvr)] 2 + 4(« + 1) sin 2 (F7r) 

^ tanh ( -3E B s ) + - ( F - -) tanh ( -/3m 

(2F-l)£ BS + m V2 J 3 V 2; V2 



Unlike L(T) and S(T), 0(T) is finite at half-integer values of e$: 
0(T) | F= i = --^ — ([1 - sgn(cos 0)] tanh (}-(3msmQ] + 



sin 20 ? A, tanh 2^^ 



2vr y - 1) v-sin^0 



(4.19) 



(4.20) 



In the cases of A = and A 1 = 0, Eq. ()4.19|) takes form 

O(T) = --Ji—sFQ. - F) (f - ± ± tanh Q/3m") , = ±a| (mod2vr) . (4.21) 

Note that relation (|4.18|) at zero temperature was first obtained in R el. 2.3.1 • and expression 
(|4.19|) at zero temperature takes form (see Ref. [22]) 



O(0) 



47r|m 



sF(l — F) 



v F A - v l - F A~ l 



00 

1 / i\ 1 r dv 

3\ ~2~) + 2^J v^/v^T v F A + 2 + v l - F A- 1 
1 



. (4.22) 



In the high-temperature limit we get 

.2 f 1 



OCT - 00) = sF(l - F) \ l - [1 - sgn(A)] 

8nk B T {2 {2F — l)E BS + m 3 



sin(F7r) /" ofat 



u A - u^A- 1 - 2mcos(Ftt) 



rr J u \u F A — v}~ F A~ l + 2 cos(-F7r)] 2 + 4(u + 1) sin (i<V) 




(4.23) 



5 Correlations 

Similar to Eq. (j3.2|) . we get the thermal correlations of fermion number with orbital angular 
momentum 

A(T;L,N) = \j ^sech 2 Q/fcA Tr A(H - u>)-\ (5.1) 
c 
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and fermion number with spin 

A(T;S,N) = ^j ^sec^Q^TrS^-o;)- 1 , (5.2) 

c 

where contour C is defined as above. Using Eqs. (|3.28|) and (|3.29|) and deforming the contour 
around the cuts and poles on the real axis, we obtain the following expressions for the 
correlations as real integrals: 



oo 

. /rr , f ~ N s sin(F7r) f du , 9 f 1 n , \ 

A(T; L, N) = A L / _ sec h 2 -pmy/u + T 

2tt J u \2 J 



X 



X 





F 2 u F A - (1 - Ffu 1 - F A- 1 - u[l - 2F(1 - F)] cos(Ftt) 
[u F A - u 1 ~ F A- 1 + 2 cos(Fvr)] 2 + A(u + 1) sin 2 (F7r) 



s r , fl - 2F(1 - F)]E BS + (2F - l)m 
8 [1 - Sgn{A)][ \2F- l)E BS + m S6Ch \^ EBS 



+ -(F--lF(l-F)sech 2 (-/3m), (5.3) 



and 



oo 

. /rr , a ,TrN s sin(i<V) /" , 9 fl _ , -\ 

A(T; 5, AT) = i '- / _ S ech 2 -pm^u + 1 x 

4vr y it V 2 / 



o 

Fu F A - (1 - F)^ 1 -^- 1 - u cos(Ftt) 
X - u^A- 1 + 2 cos(Ftt)] 2 + 4(« + 1) sin 2 (F7r) + 

+ - [1 - sgn(A)] Ebs + ( f F " 1)m sech 2 ( ±f}E B8 ) , (5.4) 
16 & V ;J (2F - l)E BS + m V 2 / 

where A is defined by Eq. (|4.5|) . and E BS is determined as a real root of Eq. (j4.6|) . 
It should be noted that, if one takes operator 

A' = A + e $ = -% x x d, (5.5) 

and defines the corresponding operator in the second-quantized theory, 

1 /v 



U = -J d 2 x[^ + ,A'^}., (5.6) 
then correlation A(T; I/, AQ is infinite. To see this, let us consider quantity (compare with 
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Eq.flHUZD; 



2- 



dtptr [A' Gr(r,<p-y, if)] 



2 sin(FTr) 



■x 



[n (uj + m) tan u u Kp{Kr)Kp{Kr') + (n + s)(u — m)e jF7r K 1 _ i ?(fi;r)K 1 _^(Kr / )] — 

oo 2 ( °° 

sin(F7r) f , / n 2 rr' r 2 + r' \ | , e y 



2F(1-F)tt 



u I dy exp 



2y 



2rr 



— y){ e$— y due-"[(l-F)fr F (u) + 

2/ 



(2/ -n(n + ^]+i S 



[^(y) - K^ F {y)\ \ + 2^e$K (^|r - r'\). 

(5.7) 



The last term in Eq. (|5.7|) diverges in the limit r' — > r, and this divergence can not be 
compensated by subtraction (as is the case for Eqs. (j3.18j) . (j3.20|) - (j3.22|) ). because Eq. (j3.27|) 
holds. Since this divergence is proportional to uj, it does not contribute to average L'(T), but 
does contribute to correlation A(T; L', N) yielding a term which is — | e$ sech 2 (|/3m) times 
infinity. Thus, if one accepts finiteness of correlations as physically plausible condition, then 
one has to favour gauge-invariant definition of orbital angular momentum, i.e. to choose 
A (l3~T7l) instead of A'lETSI). 

Contrary to the case of averages L(T) and S(T), correlations A(T; L, N) and A(T; S, N) 
are finite at half- integer values of e$: 



A{T;L,N) 



A{T;S,N) 



OO 

sin (26) f dv sech 2 (\j3my/v) 



327r J y/v — 1 t) — sin 
1 

sin , , „ s , ,o/l 



[1 — sgn(cos 0)] sech I -/3m sin© 
16 V 2 



(5i 



In the case of A = and A 1 = expressions (J5.3)) and (|5.4j) simplify: 









4) 







I 2 ( F ± - 

2 2 



1 \ 7T 

sech 2 ( 2^ m J ' ® = ~ s ~ s ~2 ( mo< ^27r), 

(5.9) 



S / „ 1 , 1 \ , 2 f 1 „ \ Ti 



and 



A(T;S,N) = --\F--±- \ sech 2 (^mj , = ± s - (mod 2vr). (5.10) 

In the limit T — > (/? — > 00) correlations ()5.3|) and (|5.4|) tend exponentially to zero for 
almost all values of with the exception of one corresponding to the zero bound state 
energy, E BS = (A = -1): 



0. 



A(0;L,JV) = --A(0;5,JV) 



2 V 2 



A^-l 
A = -1 ' 



(5.11) 
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In the high-temperature limit the correlations tend to finite values: 



s sin(Fvr) f rfw F 2 m f A - (1 - F) V"^ 1 - u[l - 2F(1 - F)} cos(Ftt) 

2^ y V [u F A - u 1 - F A~ l + 2 cos(Fvr)] 2 + 4(w + 1) sin 2 (F7r) 
o 

and 

c m s h f A „E BS + (2F - l)m , 
A(oo; 5, JV) = — [1 - sgn(A)]- 



(5.12) 



16 1 & v yj (2F- l)£ B5 + m 

oo 

ssin(Fvr) /" du Fu F A - (1 - F)u l ~ F A~ l - u cos(Ftt) 

+ 4^ y ~ - u l - F A- 1 + 2 cos(Ftt)] 2 + 4(w + 1) sin 2 (F7r) ' ^' ^ 
o 

Summing Eqs. (j5.3|) and (|5.4|) . we get the thermal correlation of two conserved observables, 
total angular momentum and fermion number, which can be recast in the form: 

A(T; M, N) = -s (f - A(T; N, N) - ^ (f - F(l - F)sech 2 Q/3mJ , (5.14) 

where 



oo 

A(T; iV, N) = f *H sech 2 ( ^m^+l) x 

2vr J u \2 J 



x 



o 

Fu F A + (1 - F)m 1_f A _1 - u(2F - 1) cos(Ftt) 
- u^A- 1 + 2 cos(Fvr)] 2 + A(u + 1) sin 2 (Fvr) 



+ - [1 - sgn(A)] sech 2 f -(3E BS \ - -F(l - F) sech 2 I -/3m 1 (5.15) 

is the quadratic fluctuation of fermion number which was first computed in Ref.j21]. Note 
that correlation A(T;M, N) vanishes at half-integer values of e$. 

Using trace identities (|3.38|) and (|3.39|) . we get the thermal correlations of total angular 
momentum with orbital angular momentum 

A(T; L, M) = -s(f- A(T; L,N) + ^[1- F(l - F)]F(1 - F) sech 2 Q/3mJ , (5.16) 

and total angular momentum with spin 

A(T; S, M) =-s(f- ^ A(T; S, N) - ±F(1 - F) sech 2 Q/3m) . (5.17) 

Using trace identity ()3.36|) . we get the thermal correlation of fermion number with induced 
flux multiplied by e 

,2 



A(T;6,N) 



27cm 



- A(T; N, N) - {F - | J A(T; S, N) - -F(l - F) sech 2 ^ ±/3m 

(5 
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or in the explicit form 



oo 

. lrT1 A ^ e 2 sin(FTr) sF(l — F) f du , 2 fl a , -\ 



x 





u F A + u l ~ F A~ l 



[u F A - u^A- 1 + 2 cos(Fvr)] 2 + 4(w + 1) sin 2 (F7r) 

^sech 2 ^ 

(5.19) 



e 2 u ,a\, sF(l-F) , 2 A OE1 A e 2 sF(l-F) 2 A 



At half-integer values of e<E> we get 

se 2 cos6 f dv sech 2 (^m^) 



A(T;0,N) 



f=\ A(An) 2 m J y/v - 1 v — sin 9 
i 

2 /1 \ „ „2 



s e 

+ 



[1 — sgn(cos O)] sech 2 ( -/5msin9 J — — sech 2 ( . (5.20) 



Correlation A(T; O, iV) vanishes in the cases of A = and A 1 = 0. In the zero-temperature 
limit we get 

A(0;O,iV) = ^ e 2 5F(l - F) A = _ x ~ (5-21) 
t 87rm 

In the high-temperature limit correlation (|5.19jl tends to a finite value 

» y » ~ n e 2 ri .... sF(l-F) e 2 sF(l-F) 

A(oc;0, N) = — [l-sgn(A)]- 



16tt l to v /j (2F - 1)£ B5 + m lQirm 

oo 

e 2 sin(F7r) sF(l — F) f du u F A + u 1 ~ F A- 1 

+ 2(2tt) 2 m ./ IT [u F A - u l ~ F A- 1 + 2 cos(Ftt)] 2 + A{u + 1) sin 2 (Fvr) ' ^' 22 > 

o 

Using trace identity (|3.4U|) . we get the thermal correlation of total angular momentum 
with induced flux multiplied by e 

A(T; 6, M) =-s(f- A(T; O, iV). (5.23) 

Thus, this correlation vanishes at half-integer values of e$. 



6 Nonnegativeness of quadratic fluctuations 

As we have seen in the previous section, a gauge-invariant definition of angular momentum 
is required by the finiteness of the correlation of angular momentum with fermion num- 
ber. In the present section we shall consider further restrictions which are imposed by the 
nonnegativeness of quadratic fluctuations. 



21 



With the use of trace identity (|3.37|) . the thermal quadratic fluctuation of total angular 
momentum, Eq. (j3.2j) . is expressed through that of fermion number: 



A(T; M, M) 



A(T; N, N) 



-F 2 (l 



FfseclT -/3m 



(6-1) 



where A(T; N, N) is given by Eq. (|5.15|) . The thermal quadratic fluctuation of fermion 
number was analyzed in detail in Ref.j21]. In particular, in the high-temperature limit 
we get 



A(oo; N,N) 



f; (1 - 








>• 









71 



^ s- (mod2?r) 



9 



7T 

s— (mod 2n) 



7[ 



^ -s- (mod27r) 



6 



7T 



0<F<- 



< F < 1 



(6.2) 



— s— (mod27r) 



which is obviously positive. Substituting Eq. (j6.2|) into Eq. (j6.1|) . we get 



A(oo; M, M) 



1(1 



-F 



la 



F — - 
2 

2 



71 



^ s- (mod27r) 







7T 



s— (mod27r) 



< F < 



71 



^ -a- (mod27r) 



Ff 



Ff 







7T 



— s— (mod27r) 



- < F < 1 
2 ~ 



fmod7r x 



(6.3) 
only, and is 



which is positive at < F < 1 - 2" 1 / 2 and 2" 1 / 2 < F < 1 ^ | 
nonpositive otherwise. 

The cause of negativeness of A(oo; M, M) is quite understandable from the mathematical 
point of view. Although operator J 2 , as a square of self-adjoint operator J, is nonnegative 
definite, appropriate spectral density Tj2(E) might be not, and the latter results in negative- 
ness of A(T; M, M), see Eq. (jl.22|) . Nonnegativeness of Tj2(E) is rooted in the procedure of 
taking a functional trace of the resolvent kernel, with regularization and renormalization in- 
volved in the case of infinite space: initially, the trace is positive but divergent (see Eq. (j3.21|) 
at r' — > r), the comparative trace which corresponds to the case of absence of the defect 
is also positive but divergent (see Eq. (j3.25J) at r' — > r), then the difference of the above 
two traces appears to be finite but not positive. Namely the same is the mechanism of the 
appearance of the negative vacuum energy density, which is widely known as the Casimir 
effect the vacuum energy density in an infinite space bounded by two parallel plates 
is positive but divergent, the vacuum energy density in an infinite unbounded space is sub- 
tracted, and the result is finite but negative. There is a physically plausible interpretation 
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of negativeness of the vacuum energy density, linking it to a force of attraction between two 
plates. Returning now to negativeness of the thermal quadratic fluctuation, we have not 
found, up to now, any physically plausible interpretation of this mathematically feasible ef- 
fect. Thus, we have to stick to the paradigma that the quadratic fluctuation of the physically 
meaningful observable is to be nonnegative. 

Let us recall that, in the case of planar rotationally symmetric system considered in the 
present paper, we have two conserved observables which are fermion number and total an- 
gular momentum. Fermion number is defined uniquely: operator N in the second-quantized 
theory corresponds to unity operator / in the first-quantized theory. As to total angular mo- 
mentum, situation is different: so far we have chosen operator M corresponding to J (|2.10|) . 
but, equally as well, the choice can be a superposition of M and N , M + EN, see Eq. (l2.16|) . 
with E being a function of the parameters of the vortex defect. We shall use this ambiguity 
and fix H by the requirement that the quadratic fluctuation of modified (improved) total 
angular momentum behave qualitatively in a similar manner as that of fermion number. 

First, it is straightforward to get general relation 



A(T; M + EN, M + EN) 



F - 



-EIF 



A(T; N,N)- 



F(l-F)sech 2 ( ^pm ) . (6.4) 



Then, one should note that, without a loss of generality, all possible values of E can be 
restricted to interval < |H| < 1/2, since shift E — > E ± 1 yields the same spectrum of 
the angular momentum operator in the first-quantized theory. Also, in view of the Bohm- 
Aharonov effect [T^], we suppose that E depends on fractional part of e$ rather than e$ 
itself, or, in other words, it depends on -F (J3.10j) . Since F ^ 1 — F as s — > — s, parameter E 
has to depend on s(F — 1/2). Taking all the above into account, we fix E = Ep, where 



1 f r, 1 

2 SSgn F "2 



1 

2' 



T 2 
2 



and define the improved total angular momentum operator as 

R = M + E F N. 
The corresponding quadratic fluctuation takes form 



(6.5) 



(6.6) 



A(T;R,R) 



F- 1 - 
2 



i) MT:X.N) + 



+ 



"1 




1 




F - 




3 




2 



-F(l-F) 



F(l - F)sech z ( ^/3m 



(6.7) 



The high-temperature limit of Eq. (|6.7J) is positive, and the behaviour of A(T; R, R) at finite 
temperatures is qualitatively the same as that of A(T; N, N). 
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7 Summary and discussion 



In the present paper we continue a study of the properties of an ideal gas of two dimensional 
relativistic massive electrons in the background of a static magnetic vortex defect, which was 
started in Ref . |24j . We find that this system at thermal equilibrium acquires, in addition 
to fermion number considered in Ref. £24] , the following nontrivial characteristics: orbital 
angular momentum (|4.3|h spin (|4.4|) . and induced magnetic flux (times e) ([4. 19)1 . The local 
features of the field strength in the interior of the vortex are exhibited by self-adjoint exten- 
sion parameter O which labels boundary conditions at the location of the vortex defect, and 
arbitrary values of vortex flux $ are permitted; our results are periodic in with period 2ir 
at fixed $ and periodic in $ with period e~ l at fixed O. Orbital angular momentum and spin 
are odd and induced flux is even under transition to the inequivalent representation of the 
Clifford algebra (s — > — s or m — > —m). In the zero-temperature limit the results of Refs.[23j 
and [22] are recovered, see Eqs. (j4.13|) . (j4.14|) and ()4.22j) . whereas in the high-temperature 
limit all averages vanish, see Eqs. (j4.15|) . (j4.16j) and (|4.23|) . 

The key point in this study is played by Section 3, where the appropriately renormal- 
ized traces of the resolvent operator are obtained, see Eqs. (|3.28|) - (|3.33|) . Thermal averages 
are then computed as integrals over a contour in the complex energy plane, see Eqs. (j4.1j) . 
()4.2|) and (j4.17|) . Moreover, the knowledge of the resolvent traces allows one to compute 
also thermal correlations of conserved and nonconserved observables and thermal quadratic 
fluctuations of conserved observables. In particular, we have computed the correlations of 
fermion number with orbital angular momentum, spin and induced flux multiplied by e, see 
Eqs. (J5.3|) . ()5.4|) and ()5.19j) . These correlations vanish at zero temperature unless the bound 
state energy in the one-particle spectrum vanishes (A = —1), see Eqs. ()5.11j) and (j5.21j) . The 
high-temperature limits of the correlations are given by Eqs. fl5.12|) . ()5.13|) and ()5.22j) . Note 
that correlations (|5.3|) and (J5.4j) are even and correlation (J5.19j) is odd under transition to 
the inequivalent representation of the Clifford algebra (s — > — s or m — > — m). 

It should be emphasized that, owing to the gauge-invariant definition of orbital angular 
momentum (j2.17|) . the correlation of orbital angular momentum with fermion number is 
finite. As it is shown in Section 5, another definition of orbital angular momentum (see 
Eq. (j5.5j0 results in the infinity of the appropriate correlation, which can be regarded as 
unphysical. 

To illustrate the behaviour of an average and a correlation as functions of the boundary 
parameter B, we depict these quantities on Figs. 1-3 for one observable, induced flux multi- 
plied by e (given by operator O (|1.14|) with Q (|2.19jl ). at several values of the vortex flux. 
Here quantities s\m\e~ 2 0(T) and sme~ 2 A(T; O, N) are along the ordinate axes, and quantity 
s07T _1 is along the abscissa axes. Values (fcgT/|m|) = 5 _1 , 1, 5 correspond to two dashed 
(with longer and shorter dashes) and one dotted lines, and values T = and T = oo corre- 
spond to solid lines; the latter cannot lead to confusion, since, as it has been already noted, 
the average at T = oo vanishes everywhere, while the correlation at T = vanishes almost 
everywhere with the exception of one point (A = —1). Our plots correspond to three values of 
F (I3.1(J|) from interval < F < 1/2, whereas interval 1/2 < F < 1 can be considered by tak- 
ing in view that s\m\e' 2 6{T) -> -s\m\e- 2 6(T) and sme~ 2 A(T; 0, N) -> sme" 2 A(T; 6, N) 
at F -> 1 - F and -> -0. 

As is seen from Figs. 1-3, the average at zero temperature is characterized by a jump with 
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a cusp at the point corresponding to the zero bound state energy (A = — 1). As temperature 
increases, this jump is smoothed out, while extremum evolves close to = s| (mod 2tt) 
(A- 1 = 0) in the case of < F < 1/2 and to 6 = -sf (mod 2tt) (A = 0) in the case 
of 1/2 < F < 1. As temperature departs from zero, the correlation develops a maximum 
at A = —1, which persists in the case of F — 1/2 and is shifted to the left in the case 
ofO < F < 1/2 and to the right in the case of 1/2 < F < 1, with further increase 
of temperature. At non-zero temperature the correlation becomes negative at cos O > 
(A > 0), i.e. in the case of absence of the bound state in the one-particle spectrum. 

Due to rotational invariance of the system considered, there is an additional to fermion 
number conserved observable — total angular momentum. It is natural to take in this 
capacity the sum of orbital angular momentum and spin: M = L + S. The appropriate 
average is related to the average fermion number, see Eqs. (j4.8|) and (J4.9|) . and the correlation 
with fermion number is related to the quadratic fluctuation of fermion number, see Eqs. fl5.14|) 
and ()5.15|) . The latter relations are a consequence of trace identity ()3.35|) . and, using trace 
identities ()3.38|) - (j3.40J) . we get the correlations of total angular momentum with orbital 
angular momentum, spin and induced flux times e, see Eqs. (j5.16|) . (j5.17|) and ()5.23|) . Using 
trace identity (|3.37j) . we get the quadratic fluctuation of total angular momentum (jbM|) . 
However, the negativeness of this fluctuation signifies that the above defined total angular 
momentum can not be regarded as a physically meaningful observable. 

To remedy the situation, in Section 6 we introduce improved total angular momentum 
as a sum of naive total angular momentum and fermion number with the dependent on the 
vortex flux coefficient, see Eqs. (|6.5|) and the quadratic fluctuation of the improved 

observable is given by Eq. (|6.7|) . On Figs. 4-6 we plot averages and quadratic fluctuations of 
naive (M) and improved (R) total angular momenta for several values of F from interval 
< F < 1/2; interval 1/2 < F < 1 can be considered by taking in view that the averages 
multiplied by s and the fluctuations are invariant under F — > 1 — F and O^—O. As 
before, two dashed (with longer and shorter dashes) and one dotted lines correspond to 
values {ksT /\m\) = 5 -1 , 1, 5, and solid lines correspond to values T = 0, oo. In the F / 1/2 
case, both fluctuations at extremely small temperatures possess a peak at A = — 1 which 
is smoothed out as temperature increases; incidentally, the minimum evolves to the left of 
point A^ 1 = in the case of < F < 1/2 and to the right of point A = in the case 
of 1/2 < F < 1. In contrast to the .^-fluctuation, the M-fluctuation is obviously negative 
in the vicinity of this minimum, see Figs. 4 and 5. In the F — 1/2 case, the behaviour 
of two fluctuation is completely different, with the M-fluctuation being independent of O 
and negative, see Fig.6 and explicit expression (|6.1|) . The behaviour of the /^-fluctuation 
is qualitatively the same as that of the fermion number fluctuation studied in detail in 
Ref . |24j . Thus, we conclude that a physically meaningful observable is the improved total 
angular momentum (R) rather than the naive one (M). 

At this point it is appropriate to discuss a very tiny effect which is common for fluctuations 
of fermion number and improved total angular momentum. Namely, both fluctuations at 
rather small temperatures become negative with extremely small absolute values in the case 
of absence of the bound state in the one-particle spectrum, i.e., at cos© > (A > 0). Fig. 7 
illustrates this fact for the -R-fluctuation. To be more specific, in the F ^ 1/2 case, this 
fluctuation at T = \m\/ (5ks) attains value — 10~ 4 at F = 0.1 (F — 0.9) and value —3 x 10~ 4 
at F = 0.3 (F = 0.7) in region (mod27r) < < sf (mod27r) (-sf (mod27r) < < 
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(mod27r)); in the F = 1/2 case, the fluctuation at T = \m\/kB attains value —27 x 10~ 4 
in region — s| (mod27r) < < (mod27r). 

Let us recall that boundary parameter B is introduced as a parameter providing a self- 
adjoint extension of the Dirac Hamiltonian in the background of a pointlike vortex defect. 
Thus the mathematical requirement of self-adjointness of the Hamiltonian is consistent with 
all values of O. However, the physical requirements may restrict the range of O. As we see, 
this really happens: nonnegativeness of the fluctuations of fermion number and improved 
total angular momentum is consistent with values of G from the range corresponding to the 
case of existence of the bound state in the one-particle spectrum, i.e. cos 6 < (A < 0). 

To conclude, we reiterate that a comprehensive study of thermal characteristics of the 
planar quantum fermionic system with a pointlike vortex defect has allowed us to reduce the 
mathematically acceptable set of boundary conditions at the location of the defect to the 
physically acceptable one. It is found that angular momentum of the system is defined as 
shifted by 1/2 from its naive value. 
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Appendix A 

The radial components of the resolvent kernel of H ()3.3|) are presented in the following way, 
type 1 (7 = s(n — n ) > 0): 



VK 



a n (r; r') = ^-(u + m) 6{r - r')H^ F (kr)J^ F (kr') + 6{r' - r)J^ F {kr)H^ F {kr') , (A.l) 



r(i) 



K(r;r') = y&; 



r - r')HW(kr) J,_ F (Ar') + 9(r' - r)J l+1 ^ F {kr)H\ L J F {kr>) , (A.2) 



(i) 



in 



c„(r;r') = —{u - m) 



6(r - r')Hil\_ F (kr)J l+1 ^ F (kr') + 6(r' - r)J l+1 „ F (kr)H^[_ F (k 



(A.3) 



d n (r;r') = ^k 



9{r - r')H^ F (kr)J l+1 ^ F (kr') + B{r' - r)J^ F {kr)H\ l ^ F {k 



type 2 (/' = -s(n - n ) > 0): 



(A.4) 



a n (r; r') = ^-(a; + m) 



r-r')HV F (kr)J v+F (kr') + e(r' - r)J v+F (kr)H\^ F (kr') , (A.5) 



(i) 



b n (r;r') 



— k 
2 



0(r-r / )^ 1) 1+F (A;r)J, +F (A;r , )+^(r / -r)J^ 1+F (A;r)^ F (A;r , ) , (A.6) 



(i) 
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lit 



c«(r; r') = ^ (u-m) [0(r - r')H^ 1+F (kr)J^ 1+F (kr') + 9(r' - r)J v . 1+F (kr)H^ 1+F (kr 



" (i. r >\ 
(A.7) 



dn(r; r') = -y/c [e(r - r')H^ F (kr)J V - 1+F (kr') + 0(r' - r)J v+Fv kr)H^_ l+F {kr')\ ; (A.8) 
type 3 (n = n ): 

' — — {o{r - r , )i/i 1 l(A;r)[sin^J_ F (A;r , ) + cos u u J F {kr')] + 

+ 0(r' - r)[sini/ w J_ F (Ax) + cos i/ w J ' F (kr)]H m F (kr')} , (A.9) 



bn (r;r') = 



Z7l 



9{r - r^H^pikr^smVuJ-Fikr') + cos v u J F (kr')}+ 



2 sin u w + cos h , LU e' lF7T 

+ 6(r' - r)[sinz/ w Ji_ F (A;r) - cosz/ w J_i +F (A;r)]iji^(A;r')} , (A.10) 



in uj — m 



j#(r — r')ifi^ F (/cr) [sin i/^Ji-i^&x') — cos v w J-i +F (kr')] + 



2 sin v w + cos ^e !Fl 

+ ^(r'-r)[sin^J 1 _ F (A;r) -cosz/ tJ J_ 1+F (A;r)]/j; i J F (A;r')} , (A.ll) 



lit 



k 



. — {0(r - r^H^kr^smu^Mkr') - cos^J_ 1+F (kr')}+ 

2 sin v w + cos ^e ,t,r I 

+ 9(r' - r)[amv u J- F (kr) + cos v u J F (kr)]H[%(kr') } . (A.12) 



Here /c = \fu 2 — w? and a physical sheet is chosen as < ArgA; < n (Irak > 0), 9(u) = 
|[l + sgn(u)], Jp(u) is the Bessel function of order p, H^\u) is the first-kind Hankel function 
of order p, and 



(A.13) 



k2F / w«, ,M-2F r ( 1 - i? ) @ 

taiii^ = — — sgn(m)(2m) tan Is- + - 

cj + m 1 \r ) \ I 4 

In the absence of the vortex defect the radial components take form: 
«n(r; r')| e$=0 = ^ + m) [0(r - r')H^\kr)J n (kr') + 9(r' - r)J n {kr)H^\kr')] , (A.14) 



bn(r;r')U =0 = ^k 8(r - r')H ( ^ s (kr)J n (kr') + 6{r' - r)J n+s (kr)H^(kr') , (A.15) 



J\ Zt(1) 



27 



27T 



c n{r; r')| e$=0 = tt(^ - m) 9{r - r')H^l s (kr)J n+s (kr') + 6(r' - r)J n+s (kr)H^ s (kr') 



(!) 



(A.16) 



d n (r; r')\ 



m . 



e<I>=0 



-k 



9{r - r')H^(kr)J n+s (kr') + 0(r' - r) J n (kr)HSl a (kr')\ . (A.17) 



Note that all radial components behave asymptotically at large distances as outgoing 
waves. 

Appendix B 

Let us consider quantity (j.lllj) . The contribution of the regular (types 1 and 2) components 
of the resolvent kernel is given by expression 

[{n — e$)a„(r; r') + (n + s — e$)c„(r; r')\ = 

= s [{l + m)I^ F (Kr)K^ F (Kr') + (Z + 1 - F){u - m)I l+1 „ F (Kr)K l+1 - F (Kr')} - 

i>i 

-« K Z ' + + WM^'+iK™"') + (/' - 1 + F)(u - m)I v _ l+F (Kr)K v _ l+F (Kr')} , 

(B.l) 



V>1 



where k = —ik and we used relations 



J p (iKr) = e ip ^I p (Kr) and H {l \iKr') = —e~ ip ^K i 

ITT 



\Kr) 



which are valid at Ren > 0; here I p (u) is the modified Bessel function of order p, and 



7T 



K p( U ) = o • / \ L--/' 

2 sm(p7r) 



[/_» - /»] . 



Using relation (see, e.g., Ref . [5Dj ) 

oo 

I p {Kr)K p {Kr ) = - J — exp 



2y 



2 + r' 2 \ 



^(Z + = [/ p (y) + I P+ i(y)] 



(B.2) 
(B.3) 



z>i 



we perform summation in Eq. (|B.1J) and get in the case of Ren > \Im n\: 

[(n — e$)a n (r; r') + (n + s — e$)c n (r; r')] = 



s sin(F7r) 



7T 



uj I dy exp 



2?7 



2rr' 



— y [AT F (y)-*i-F(y)] + 



+ sF(cj + m)I F (Kr)K F (Kr') - s(l - - rr^h^K^Kx^Kr') . (B.4) 
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The contribution of the irregular (type 3) components is given by expression 

/ t\ / a / t \ / /\ 2ssin(F7r) 

(n - e$)a no (r; r ) + (n + s - e$)c no (r ; r ) = — x 

7r(tan i/ w + e* r "j 

x [-F(u; + m) tani/ w # i? (/cr)tf Ji .(«r / ) + (1 - F)(u - m)e iFw K l _ F (Kr)K 1 _ F (Kr')] - 

- sF(u + m)I F {Kr)K F {Kr') + s(l - - m)h- F {Kr)K X - F {Kr') . (B.5) 

Summing Eqs.flBl} and flB.5H . we get Eq. ffTTTjl . 

Computation of other quantities, Eqs. (|H.12j) - (|H.lfjjl . is similar to the above. As an illus- 
tration, let us scrutinize the procedure of computation of the last quantity, Eq. (j3.16j) . which 
is the most tedious one. 

The contribution of the regular components is 



e 2 

— sr 

47T 

n^no 



(n-e® + ^ [b n (r; r') + d n {r- r')] = 

Kr fe ( Z " F + I) [Ii+i-F(^)K^ F (Kr') - I^ F {Kr)K l+1 ^ F {Kr')} 
-Y,(l' + F-^) [Iv-i + f{kt)K v 



i. 

477 ' 



: v+F (Kr') - J r+i r(/«r)^/_ 1+F («;r / )] > . (B.6) 



Using recurrency relation 

ud u I p (u) = ±pl p {u) + ul p ±i(u) 



we get 

.,2 



e 

— sr 

47T 

n^no 

e 2 



5^ (n-e$+ ^sj [b n (r; r') + d n {r; r')] = 
^{(l-F + l) [(/- F) ^M^M + (/ + 1 - F)I l+1 ^ F {^)K l+1 ^ F {Kr')} + 
+ (l + F - [(I + F)I l+F {Kr)K l+F {Kr') + (/ - 1 + F)/ w+F («r)^_ 1+F (V)] } - 

- + F - [^ +F (/tr)ir /+F («:r / ) - J,_ 1+F (Kr)tf,_i +F («r')] } = 
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16tt 



CO 

d v ex P 

i>i i 



2y 



2 _|_ \ r 3 i i 



3 13 3 1 

+ ^h+2-F(y) + 2 J '-i+H^) + ^l+x+Fiy) + 2 Il - 2+F ^ + 2 Il+F ^ + 

+ (1-1- F)k^ F (y) - (I + 1 - F)I l+1 _ F {y) + (I - F)/,_ F (j/) - (/ + 2 - F)/ i+2 _ F (y) + 

+ (/-! + F)Ii-i +F (y) -{1 + 1 + F)I l+1+F (y) + (/ - 2 + F)I^ 2+F {y) - (/ + F)I l+F {v) + 



+ 



f^^'f'Y'^ <y^ ry^^ 1 



+ 



(1 - F)h_ F {y) + FI F (y) + -/^(y) + 



+ -^h+Fiy) + Ii- 1+F {y) 



, (B.7) 



where the second equality is obtained with the use of representation ()B.2jl and recurrency 
relation 

2pl p (u) = u [I p -x{u) - I p+ x{u)\ . 

Now the summation over I can be performed with the use of Eq. (jB.3|) and relation (see 
Ref.gUj) 



1=1 



p 



due u I p {u) - yl p {y) - yl p+ i{y) 



L o 



Rep> -1. 



(B.8) 



Using again the recurrency relation and the relation between the Macdonald and the modified 
Bessel functions, we get 



e 

An 



sr ^ (n - e$ + j [b n (r; r') + d n (r; r')] 
e 2 sin(F7r) 



8F(1 - F)vr 2 
y 



dy 1 + 



cxp 



2y 



2rr 



— V x 



x | e y J due- u [(l - F)K F (u) + FK^ F {u)] + (2F - l)y[K F (y) - K x _ F {y)\ | + 

2 sin(Fvr) 



•57T 



exp 



2y 



2rr 



—y 



71 



[{l-F)K F {y)+FK^ F {y)\- 



- [F - - [FI F (y) - (1 - F)h_ F (y)} + 77 



1 / /t 2 rr' r 2 — r' 2 



^ - 2 ) [If(v) ~ h-F 



(B.9) 
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The contribution of the irregular components is 

e ' ( * 1 \ r, , a J , /m e2 ( F - I) sin ( F?r ) 

x [tan^is: 1 _ F (/tr)K" F (/tr / ) -e^ J f^ F (/^r)i^ 1 _ F (/^r , )] + 

+ £ ( F ~ Kf it-i+F^KpiKr') - I„ F ( K r)K^ F ( K r')}. (B.10) 

The terms containing the 7 p -functions are cancelled in the sum of Eqs. (|B.9)l and (jB.lOJ) . 

y oo oo 

Decomposing the integral over u as J = J — J , we get the following expression for this sum: 

y 



o OO 



—sr 

4:71 

n=— oo 



^ (n - e$ + j [b n (r;r') + d n (r;r')} 



°° 2 

e 2 f , f K 2 rr' r 2 — r' \ / K 2 rr' (r — r') 2 

~ 4^ J y \ + 4y 2 + Arr' ) 6XP \ 2y 2rr' V 

o 

oo 

e 2 sin(F7r) /" / K 2 rr' r 2 — r' 2 \ / /? 2 rr' r 2 + r' 2 \ 

" ml - F)n 2 J dy \ + 6XP \~y V ) * 

o 

x^jdu e" u [(l - F)^(w) + FKi- F {u)} - (2F - l)y[K F (y) - K x . F {y)\ j + 

°° 2 

e 2 sin(F7r) /" , / K 2 rr' r 2 + r' \ r/1 T ^ , N _, T ^ , Nl 
+ 4? ; 2 y dy exp — [(l - F)K F {y) + F K x . F {y)\ - 

o 

e 2 (F - |) sin(FTr) 



2vr 2 (tan^ + e iFw ) 
Using relation (see, e.g., Ref . [3T] ) 



[tfwv u K 1 - F (Kr)K F (Kr') - e iF7T K F {nr)K^ F {nr')] . (B.ll) 



oo £ 

y dyy'- l exp{-py - gy" 1 ) = 2 Q * K s {2^/pq), 
o 

we express the first integral over y in Eq. (|B.llJ) through the Macdonald function and get 
En. (|5^2jl . 

In the case of the absence of the vortex defect one uses summation formulae 

oo oo 
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Thus, in particular, we get 



e 2 °° / 1 \ 

— sr ^2 ( n + 2 s ) M r ; r ') + dn(r;r')] = 

n=-oo V ' e$=0 

6 2 00 / j\ 

= "ivr^ S ( S ™ + 2 ) [Wl( Kr )^n(^0 - / S n(^)^ S n+l(^')] 



n=— 00 

00 



e 2 00 / j \ 
— ^ ( + - j {sn/ sn (Kr)K sn (/tr') + (sra + l)/ sn+1 (/«r)i ; C sn+ i(Kr / ) 

n=— 00 ^ ' 

- rd r [I sn (Kr)K sn (Kr') - I sn+1 (Kr)K sn+1 (nr')]} = 

OO 

e 2 ^ y / K 2 rr , r 2 + r /2 x / x x 



n=— 00 



H J [ J sn(^/) ~ Wl(?/)]} = 

e 2 r / « 2 rr' r 2 — r' \ / K 2 rr' (r — r') 2 \ 

= syH + ^ + ~^)^r~*T~ a^V' (B12) 





which is Eq.(JS2SJ). 
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